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We consider a space-homogeneous gas of inelastic hard spheres, with a 
p I ■ diffusive term representing a random background forcing (in the framework 

of so-called constant normal restitution coefficients a £ [0, 1] for the inelas- 
f-n | ticity). In the physical regime of a small inelasticity (that is a £ [a*, 1) for 

some constructive a* S [0, 1)) we prove uniqueness of the stationary solution 
for given values of the restitution coefficient a 6 [a*,l), the mass and the 
momentum, and we give various results on the linear stability and nonlinear 
stability of this stationary solution. 
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1 Introduction 



1.1 The model and main result 



We consider the steady states of the spatially homogeneous inelastic Boltzmann 
equation for hard spheres with thermal bath forcing. More precisely, we consider a 
gas which is described by the distribution density of particles f — f(v) > with 
velocity v G M. N (N > 2) and such that / satisfies the stationary equation 



(1.1) Q a (/,/)+rA w / = in 
with a constant r > 0, and given mass and vanishing momentum: 

(1.2) / fdv = pe(0,oo), [ fvdv = 0, 



The term fiA v f, with constant \i G (0, oo), represents the effect of the heat bath. 

The quadratic collision operator Q a (f,f) models the interaction of particles by 
means of inelastic binary collisions with a constant normal restitution coefficient 
a G [0, 1) (see [6j [lj [U [9]) which preserves mass and momentum but dissipates 
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kinetic energy. We define the collision operator by its action on test functions, or 
observables. Taking ip — ip(v) to be a suitably regular test function, we introduce 
the following weak formulation of the collision operator 



(1.3) / Q a (g, f)ipdv = / / / b\u\g*f - da dv dv*, 

Jr n J J JlR^xIR^xS^- 1 

where we use the shorthand notations / := f(v), g* := g(v*), ip' := ip(v'), etc. Here 
and below u — v — v* denotes the relative velocity and v', v[ denotes the possible 
post-collisional velocities (which encapsule the inelasticity of the collision operator 
in terms of a). They are defined by 

/-, „\ . w v! . w u' 

^ v = 2 + r v * = ^~r 



with 



1 — a\ ( 1 + a 
w = v + v*, u' — ( — - — 1 u + I — - — ) \u\ a. 



We also introduce the notation x = x/\x\ for any x € R , 0. The function 
b = b{u ■ a) in ( 11. 31) is (up to a multiplicative factor) the differential collisional 
cross-section. We assume that 

(1.5) b is Lipschitz, non-decreasing and convex on (—1, 1) 
and that 

(1.6) 3b m ,b M e (0,oo) s.t. Vac e [-1,1], b m < b(x) < b M . 
Note that the "physical" cross-section for hard spheres is given by (see [Hill]) 

N — 3 

(1.7) b(x)=b' (l-x)-—, 6' e(0,oo), 



so that it fulfills the above hypothesis (jl.5|1.6p when N = 3. These hypothesis are 
needed in the proof of moments estimates (see [8j Proposition 3.2] and [91 Proposi- 
tion 3.1]). 

We also define the symmetrized (or polar form of the) bilinear collisional operator 
Q a by setting 

Q a (g, h) ip dv = - III b\u\ g*h da dv dv*, 



{ JRN 2 J J J R N X RN xS N-l 

with A^, = (ip' + ip* -ip -ip*). 

In other words, Q a {g, h) = {Qa(g, h) + Q a (h, g))/2. The formula (11.31) suggests the 
natural splitting Q a = — Q~ between gain and loss part. The loss part Q~ can 
be defined in strong form noticing that 




(Qa(gJ): ^) = I I I b\u\g*fi/)dadvdv* =: (f L(g), ip), 
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where (■, ■) is the usual scalar product in L 2 and L is the convolution operator 



(1.9) L{g){v) — (fo | ■ | *g)(v) — b / g{v*) \v — v*\ dv*, with b = / b(ai) da. 

Jrn Js n - 1 

In particular note that L and Q~ = Q~ are indeed independent of the normal 
restitution coefficient a. 

As explained in [6], the operator (11.31) preserves mass and momentum, and since 
the Laplacian also does so, the mass and momentum of a stationary solution can 
be prescribed. However energy is not preserved neither by the collisional operator 
(which tends to cool down the gas) nor by the diffusive operator (which harms it 
up). Competition between these two effects gives rise to a balance equation (obtained 
after having multiplied equation (11.11) by \v \ 2 and integrated) which reads as follows 

(1.10) (l-a 2 )D £ (f) =n2Np. 
The energy dissipation functional is given by 

D £ (f) -=bi f /* \u\ 3 dv dv*, 



where b\ is (up to a multiplicative factor) the angular momentum defined by 
(1.11) h := - I (1 - (u ■ a)) b{u ■ a) da. 



O JgN-1 

In order to establish (jl.lOp we have used (jl.8p and the elementary computation 

\\*(v,v*,a) = (!-(«• o-)) \u\ 2 . 

Existence and qualitative properties of the steady solutions as well of the solu- 
tions to the associated evolution equation (Cauchy theory) was done in [6] (see also 
PP). Among others, it is proved the following in these papers: 

Theorem 1.1 ([6, Theorem 5.2, Lemma 7.2], JT\ Theorem 1]) For any given in- 
elastic coefficient a G (0, 1) and any given mass p G (0, oo) there exists at least one 
solution F G «S(R ) to the stationary equation U.1\) - I[J~2]) with mass p and vanishing 
momentum which furthermore satisfies 

(1.12) F(v) > a ie - a2H3/2 Vv G R N and [ F(v) e as ^ m dv < oo, 

Jr n 

for some constants ai, a 2 , G (0, oo). . 

Our main result gives a partial answer concerning uniqueness of solutions ob- 
tained by Theorem ll.lt 

Theorem 1.2 There is some constructive a* G (0, 1) such that for any a G [a*, 1], 
and any given mass p G (0, oo) the stationary equation U.l\) - I[l~2]) admits a unique 
solution F with mass p and vanishing momentum in the class of functions given by 
the existence Theorem \l.l\ 
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1.2 Rescaled variables and elastic limit a — ► 1 

Let us start with a remark. For any solution F to the equation (11. ip and any 
constant A £ (0, oo), the rescaled function g associated to / by the relation 

g(v) = \ N f(Xv). 

satisfies 

(1.13) Q a (g,g) + X- 3 rA v g = 0. 

In order to prove (11.131) we have just use the homogeneity properties Q a (g, g)(v) = 
X^QMJ)^) andA v g = \ N+2 (A v f)(\v). 

This elementary remark shows that we may choose r £ (0, oo) arbitrarily in 
equation (11.11) and we now make the choice 

(1.14) T = r a = p(l- a), 

and denote by F a a solution to equation (II. ip with r = r a : 

(1.15) Q a (F a ,F a )+T a AF a = in R N . 

At a formal level, it is immediate that with this choice of scaling, in the elastic limit 
a — > 1, the equations (ll.ip - (ll.2l) become 

( Q 1 (F 1 ,F 1 )=0 in R N , 

(1.16) < /■ /• 

/ F x dv = p, I F lV dv = 0, < Fx £ <S(R ). 

Moreover, from (11.101) . one gets 

(1.17) 2iV(l-a)p 2 - (1 -a 2 )D e (F a ) = 0. 

Dividing the above equation by (1 — a) and passing to the limit a — > 1, one obtains 

(1.18) Np 2 -D E {F 1 ) = 0. 

It is straightforward that the only function satisfying the constraints (I1.16P and 
(I1.18P is the Maxwellian function 

(1-19) Fi := M Sl = M p Sl 

where, for any p, 9 > 0, u £ H.^, the function M P)Uj g denotes the Maxwellian with 
mass p, momentum u and temperature 9 given by 

O \v-u\ 2 

(i- 20 ) M p>u ,e(v):= j^y^e-^, 

and where the temperature 9\ £ (0, oo) is given by (we recall that b± is defined 
in (HHD) 

(1.21) 0i = o^7T / M w (^)|^| 3 ^' 



2 6f 

For instance in dimension iV = 3 we obtain ^ = (3 2 7r) 1//3 /(2 10 6 2 ) 1//3 . Moreover, in 
the particular case of the hard-spheres cross-section (I1.7P in dimension 3, we find 
b-y = 6' (47r)/3 and therefore X = 3 4 / 3 /(2 14 vrfe^ 2 ) 1 / 3 . 
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1.3 Physical and mathematical motivation 

In this short paper we shall not review the physical and mathematical on the ki- 
netic theory for granular gases. Let us only refer to some key references where the 
reader can find all the desired details: for a detailed physical introduction we refer 
to [2111]; arid for a short mathematical introduction see [14J; finally apart from the 
many works on "pseudo-maxwell molecules" (see the bibliography of [H]) let us 
only mention the works directly related to inelastic hard spheres (freely cooling or 
in thermal bath): [61 [Tj, [U [9], [10] (in particular see the introductions of these papers 
for some mathematical discussions about the models and their connection to phy- 
ics). This paper is inspired from the models and questions arised in [6] and the new 
mathematical tools developed in [TO] . 



1.4 Notation 



Throughout the paper we shall use the notation (-) = a/1 + |- | 2 . We denote, for 
any p G [1, +oo], q G M and weight function uj : WL N — > R + , the weighted Lebesgue 
space L p {oj) by 



If q (u) := {/ : R N 
with, for p < +oo, 



measurable ; ||/||zg( w ) < +00} 



\f{v)\ p (v) pq Lu(v) dv 



1/p 



and, for p = +00, 



L»(R") = SU P \f( v )\ (v} q uj(v). 



We shall in particular use the exponential weight functions 

(1.22) m = m s>a (v) := e~ aW for a G (0, 00), s G (0, 1), 

or a smooth version m(v) := e - ^'''' 2 ^ with ( G C°° is a positive function such that 
£(r) = r s l 2 for any r > 1, with s G (0, 1). 

In the same way, the weighted Sobolev space W£' p (u) (k G N) is defined by the 
norm 

i/p 



\s\<k 



\d s f(v)\\ P L , ( 



and as usual in the case p = 2 we denote H^(uj) = Wg ,2 (ui). The weight ui shall be 
omitted when it is 1. Finally, for / G L\ k) with k > 0, we introduce the following 
notation for the homogeneous moment of order 2k 



mfc(/) := 



/ \v\ 2k dv 
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and we also denote by p(f) = m (/) the mass of /, S(f) = m^/) the energy of / 
and by 9(f) = S(f) / (p(f) N) the temperature associated to / (when the distribution 
/ has mean). For any p,S G (0, oo), u G M. N we then introduce the subsets of L l 
of functions of given mass, mean velocity and energy 

C p ,u '■= \h E L\\ I hdv = p, / hvdv = pu>, 
I Jr n Jr n J 

C P ,u,s '■= \ h G L\] I hdv = p, / hvdv = pu 1 I h \v | 2 dv = 8 > . 

I Jrn ' Jrn J r n J 

For any (smooth version of) exponential weight function m and any k, q > we 
introduce the Banach spaces 

(1.23) L^m -1 ) = L^m- 1 ) n C , and Wj'^m -1 ) = IVj'^m- 1 ) n C , - 
1.5 Summary of the results 

Our results, that we state now, deal with the evolution equation 



;i-24) ^ = Q a (f, f) + r a A/, /(0, •) = / in G C„ l0 , 



where r a is defined in f ll.141) . and with the associated stationary equation (11. 1 5 p 
with F G Cpfi- 

Theorem 1.3 There is some constructive a* G (0, 1) such that for a G [a*, 1], and 
any given mass p G (0, oo), we have: 

(i) For any r > 0, the equation / II. lj) admits a unique non-negative stationary so- 
lution with mass p and vanishing momentum. We denote by F a the stationary 
solution obtained by fixing r = r a (defined by 

(ii) Let define F% = M p § Q 1 the Maxwellian distribution with mass p, momentum 
and "diffusive thermodynamical temperature" 9\ defined in M.21\) . The 
path of stationary solutions a — > F a parametrized by the normal restitution 
coefficient is C 1 from [a*,l] into W k ' 1 H L 1 (e a '"') for any k G N and some 
a G (0, oo). 

(Hi) For any a G [a*, 1], the linearized collision operator 
(1.25) C a h:=2Q a {F a ,h) + T a Ah 

is well-defined and closed on L 1 (m _1 ) for any exponential weight function m 
with exponent s G (0, 1) (defined in U. 22) ). Its spectrum decomposes between 
a part which lies in the half-plane {Re£ < p} for some constructive p < 0, 
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and some remaining discrete eigenvalue \x a . This eigenvalue is real negative 
and satisfies 

3 

(1.26) \i a = — ■=- p (1 — a) + Oil — a) 2 when a — > 1. 

0\ 

The associated eigenspace is of dimension 1 and then denoting by <p a = 4> a (v) 
the unique associated eigenfunction such that \\(f) a \\Ll = 1 an d 0<*(O) < 0, there 
holds (j) a G S(M N ) (with bounds of regularity independent of a) and 

(1.27) := c (M 2 -iV^i)Fi as a -> 1, 

where cq is the positive constant such that = 1. Finally one has con- 

structive decay estimates on the semigroup associated to this spectral decompo- 
sition in this Banach space (see the key Theorem \5.1\ and the following point). 

(iv) The stationary solution F a is globally attractive on bounded subsets of L\ under 
some smallness condition on the inelasticity in the following sense. For any 
p,S Q ,M G (0, oo) there exists a** G (a*, 1), C* G (0, oo) and r\ G (0, 1), such 
that for any initial datum satisfying 



0< fm G L\ nC pfi , £o , \\f in \\ L i < M , 
the solution g to jl.2J$ satisfies 
(1.28) \\ft-F a \\ L \ <e^^. 

(v) Moreover, under smoothness condition on the initial datum one may prove a 
more precise asymptotic decomposition, and construct a Liapunov functional 
for the equation {1.21$ . More precisely, there exists k* G N and, for any 
exponential weight m as defined in U.22\) and any p,£ ,M G (0, oo), there 
exists a** G (a*, 1) and a constructive functional 7i : H k * flL 1 (m~ 1 ) — > K. such 
that, first, for any initial datum < f- ui G H * R L 1 (m _1 ) R Cp,o,£ satisfying 

the solution f to {1.2J$ satisfies 

(1-29) f{t, -) = F a + c a (t) <p a + r a {t, ■), 

with c a (t) G M and r a (t, •) G Ll(W N ) such that 
(1.30) \c a (t)\ < C^ 1 , \K(t,-)\\ L i < ae (3/2) "°'. 

And second when the initial datum satisfies additionally 

the solution satisfies also 

1 1— > Ti(g(t, •)) strictly decreasing 
(up to reach the stationary state F a ). 
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Remark 1.4 All the constants appearing in this theorem are contructive, which 
means that they can be made explicit, and in particular that the proof does not use 
any compactness argument. Unless otherwise mentioned, these constants will depend 
on b, on the dimension N , and on some bounds on the initial datum but never on 
the inelasticity parameter a G (0, 1]. All the other remarks made in fTU[ Section 1.6] 
also apply here. 

1.6 Method of proof and plan of the paper 

The general structure of the proof is inspired from [10]. We have packed into the 
appendix the technical results from this previous paper which are used here also, 
as well as a technical result of lower bound on the diffusive inelastic Boltzmann 
equation. Therefore the overall structure of the proof is likely to be more visible 
than in [10J. In the whole, we mainly prove here the few points which differs from 
[TO] (due to the replacement of the anti-drift term by a diffusive term) and refer to 
[TO] for more details. 

The first main idea of our method is to consider the rescaled equations (11.151) - 
(I1.14p with an inelasticity dependent diffusion coefficient r a which exactly "compen- 
sates" the loss of elasticity of the collision operator (in the sense that it compensates 
its loss of kinetic energy). This scaling allows to prove uniform bounds according 
to a for the family of stationary solutions F a to the equation (11.151) (recall that in 
this scaling, the diffusion is evanescent in the elastic limit). 

The second main idea consists in decoupling the variations along the "energy 
direction" and its "orthogonal direction". This decoupling makes it possible to 
identify the limit of different objects as a — > 1 (among them the limit of F a ). 

The third main idea is to use systematically the knowledges on the elastic limit 
problem, once it has been identified thanks to the previous arguments. In particular 
we use the spectral study of the linearized problem and the dissipation entropy- 
entropy inequality for the elastic problem. This allows to argue by perturbative 
method. Let us emphasize that this perturbation is singular in the classical sense 
because of the addition of a (limit vanishing) second-order derivative operator, but 
also because of the gain of one more conservative quantity at the limit. 

In Section 2, we use the regularity properties of the collision operator in order 
to establish on the one hand that the family (F a ) is bounded in H°° D L 1 (m~ 1 ) 
uniformly according to the inelastic parameter a (the key argument being the use 
of the entropy functional which provides uniform lower bound on the energy of F a ) 
and on the other hand that the difference of two stationary solutions in any strong 
norm may be bounded by the difference of these ones in weak norm (the key idea is a 
bootstrap argument). This last point shall allow to deal with the loss of derivatives 
and weights in the operator norms used in the sequel of the paper. 

In Section 3, we prove that F a — > F\ when a — > 1 with explicit "Holder" rate. 
The cornerstone of the proof is again the decoupling of the variation F a — F\ between 
the "energy direction" and its "orthogonal direction" . 
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Finally in Section 4, we prove uniqueness of the profile F a for small inelasticity 
by a variation around the implicit function theorem, in Section 5 we prove results on 
the localization of the spectrum of the linearized equation for small inelasticity, and 
in Section 6 we prove some (semi)-global stability results by combining the previous 
linearized study with entropy production estimates. 

2 Estimates on the steady states 

In this section we prove various regularity and decay estimates on the stationary 
solutions (or the differences of stationary solutions), uniform as a — > 1, which shall 
be useful in the sequel. 

2.1 Uniform estimates on the steady states 

For any a G (0, 1) we consider T a the (not empty) set of all the solutions given 
by Theorem 1 1.1 1 of the inelastic diffusive stationary Boltzmann equation (11.151) with 
inelasticity coefficient a, given mass p G (0, +oo) and finite energy. More precisely, 
we define T a as the following set of functions 



We show that for any stationary solution F a G T the decay estimates, the 
pointwise lower bound and the regularity estimates can be made uniform according 
to the inelasticity coefficient a G [ao ; l)- Let us emphasize once again that the 
choice of the rescaling parameter r a = p(l — a) in (11.151) is fundamental in order 
to get uniformity of these bounds in the limit a — > 1. Let us also mention that our 
choice of scaling for the equation (11.151) is mass invariant, that is F with density 
p(F) satisfies the equation if and only if F/p(F) satisfies the equation with p = 1. 
Therefore all the estimates on the profiles are homogeneous in terms of the density 
P- 

Proposition 2.1 Let us fix «o G (0, 1). There exists a±, a 2 , a 3 , a± G (0, oo) and, for 
any fceN, there exists Ck G (0, oo) such that 

(2.1) Va G [a , 1), V F a G T a , \\F a \\ Ll ^ ai \v\^ < a 2 , \\F a \\ H k^ R Nj < C k . 

Proof of Proposition \2.1\ We split the proof into several steps. We fix a G [a , 1) 
and F a G T a for which we will establish the announced bounds. Let emphasize that 
thanks to the a priori bounds satisfied by F a all the computations we will perform 




satisfying 




For some fixed a G (0, 1), we also define 
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are rigorously justified. From now we omit the subscript "a" when no confusion is 
possible. 

Step 1. Upper bound on the energy using the energy dissipation term. We prove 
that 

(2.2) Vote (0,1] 8<p( — \ . 

From equation (I1.17P on the energy of the profile G there holds 



(2.3) (l + a)h / F F* \u\ dv dv* =2N p . 

Jr n Jr n 

From Jensen's inequality 

/ \u\ 3 F*dv* > p\v\ 3 , 
Jr n 

and Holder's inequality 

/ \v\ 3 F dv > p- 1/2 ( [ \v\ 2 Fdv) , 
Jr n \Jr n J 

we get 

{l + a)b l p 1 / 2 £ 3 ' 2 <2Np 2 
from which the bound (12.21) follows. 

Step 2. Lower bound on the energy using the entropy. We prove 

fa 2 N 2 \ 2/3 

(2.4) V«6(0,l] £>p^—-\ W ith 6 a :=||6|| L i. 

Remark 2.2 The choice of scaling we have made for the rescaled equation M.15\) 
becomes clear from this computation: it is chosen such that the energy of the station- 
ary solution does not blow up nor vanishes for a — > 1. The restriction a 6 [oq, 1), 
ao > 0, is then made in order to get a uniform estimate from below on the energy. 

By integrating the equation satisfied by F against log F we find 
Q(F, F) log F dv + p (I -a) / \ogF A v Fdv = 0. 



Then we write the first term as in [HI Section 1.4] to find 

1 r r r f w w ww \ 

-J\\ FiU log— f -— f + 1) Bdvdv.da 



1 

+ 2 




xS N ~ 



{F'Fl — FF*) B dv dv* da — p(l — a) [ dv = 

-i Jrn F 
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Recalling that x — log x — 1 > for any x > and making the change of variables 
(v,v*) in the second term we obtain 



(*',<) 



(2.5) p(l-a) 



FF„ Iwl dvdv*. 



On the one hand, from Cauchy-Schwarz's inequality 



F F* Ifl (if (if* < 



(2.6) < 



F F* (if (if* 



1/2 



F FAu \ 2 dv (if* 



1/2 



72 o^S 1 ' 2 . 



On the other hand, we compute 



< 



2VVF + P -^VF 



£ 



dv 



[ \a\vVf\ 

Jr n \ 



o pNv „ / pN 



If | 2 F du 



(2.7) < 2 



V^I 2 d „_ 2 M)! + ^l f = 2 



f e \ £ 

Gathering ijESjl . (CTl and (EH) we get 



2£ ~ 2 a 



from which we deduce (|2 



S'iew 5. Upper bound on exponential moments. There exists A, C > such that 



Va G [0,1) 



F(v)e AH dv < C p. 



We refer to [U Theorem 1] where this bound is obtained as an immediate con- 
sequence of the following sharp moment estimates: there exists an explicit X > 
such that 



(2.8) 



Vae[0,l), m k = [ G\v\ k dv <T(k + 1/2) X k/2 p. 



It is worth noticing that in [T] the Povzner inequality used in order to get ( 12. 8p is 
uniform in terms of the normal restitution coefficient a G [0, 1] appearing in the 
collisional operator and in terms of the coefficient p in front of the thermal bath 
term. Moreover the factor p comes from our choice of the scaling variables (in which 
p is involved). 
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Step 4- Uniform upper bound on the L 2 norm, and H k norms, k > 0. On the 
basis of the uniform bounds from below and above on the energy, the proof can be 
done exactly as in [TUJ Proposition 2.1], since the diffusion term plays no role in the 
energy estimates (it only helps). □ 

After the estimates from above on the profile, let us now state a pointwise bound 
from below (it is a straightforward consequence of Proposition ID . II in the appendix). 

Proposition 2.3 Let us fix «o G (0, 1). There exists as, a 4 G (0, oo) such that 

(2.9) V«6[q ,1), VF q G^ q , F a >a 3 e- a ^\ 

2.2 Estimates on the difference of two stationary solutions 

In this subsection we take advantage of the mixing effects of the collision operator 
in order to show that the L 1 norm of their difference of two stationary solutions 
(corresponding to the same inelasticity coefficient) indeed controls the if fe nL 1 (m _1 ) 
norm of their difference for any fceN and for some exponential weight function m, 
uniformly in terms of a G [oto, 1). 

Proposition 2.4 For any k > 0, there is m = exp(— a \v\), a G (0, oo) and Ck > 
such that for any a G [«o, 1) and any F a , H a G T a there holds 

(2.10) \\H a — -Fall/ffcnLUm- 1 ) < Cfc \\H a ~~ -^alU 1 - 

Proof of Proposition \2.4[ We proceed in three steps. It is worth mentioning that all 
the constants in the proof are uniform in terms of the normal restitution coefficient 
a G [«o, 1), as they only depend on the uniform bounds of Proposition 12. II and some 
uniform bounds on the collision kernel. 

The proof is based on the following three steps, which can all be proved following 
exactly the arguments in [TD"1 Proposition 2.7]. 

Step 1. Control of the L 1 moments. We prove first that there exists A,C G (0, oo) 
such that 

Vae[a ,l), / \H a - F a \e Alvl dv <C [ \H a -F a \dv. 
Jr n Jr n 

The proof can be done exactly as in the Step 1 of [TUJ Proposition 2.7]: the only 
change is the replacement of the vanishing anti-drift term by the term viscosity term, 
which both in any case are negligible in the moments estimates (note also that this 
estimate on the diffusive Boltzmann equation is proved within the paper [I], with 
constant independent on the viscosity coefficient). 

Step 2. Control of the L 2 norms. It is done exactly as in Step 2 of [TO], Proposi- 
tion 2.7], using the regularity theory of the gain term: the only difference is that 
instead anti-drift term which vanishes in the estimates, one has now a diffusion term 
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which yields good damping negative in the L 2 energy estimate (which can therefore 
be dropped). 

Step 3. Control of the H k norms. The proof can be done exactly as in Step 3 of [10, 
Proposition 2.7]. The anti-drift is replaced by the diffusion term which is also a 
good negative term in the H k energy estimates. □ 

3 Quantification of the elastic limit a — > 1 

We have the following estimate on the distance between F a and _F\ for any stationary 
solution F a . 

Proposition 3.1 For any e > there exists C e (independent of the mass p) such 
that 

(3.1) V«G[« ,1) sup \\F a -F x \\ Ll <C £ p{l-a)^ 
where we recall that F\ is the Maxwellian function defined by U.19\ )- fl.21\) . 

Proof of Proposition UTTi On the one hand, for any inelasticity coefficient a G [«o, 1) 
and profile F a , there holds from (12.51) together with Corollary IA.4I and the uniform 
estimates of Proposition 12.11 

(3.2) D HA {F a ) < D H<a (F Q ) + p 2 0(1 - a) < p 2 0(1 - a). 

On the other hand, introducing the Maxwellian function Mg with the same mass, 
momentum and temperature as F a , that is Mg given by (11.201) with u = and 
9 = £(F a )/p, and gathering (13.21) . (IA.9I) . (1A.8I) with the uniform estimates of Propo- 



sition 12.11 and interpolation inequality, we obtain that for any q, e > there exists 



C q)£ such that 



(3.3) Vote [«„,!) \\F a - M e \\ 2 L r <C q , £ p 2+e (l-a). 
Next, from (12.31) . we have 

b\ [ [ F a F a *\u\ 3 dv dv*—2 p N f F a dv = (l—a)— f f F a F a * \u\ 3 dv dv* 

JrnJ-^n jRiv 2 JrjvJkjv 

and then 

(3.4) |tf (0)| < d \\F a - M e \\ L i + C 2 p 2 (1 - a), 

where we have used that F a and Mg are bounded thanks to Proposition 12.11 and we 
have defined 



(3.5) 



m(6)=2pN I Mgdv-h / MgMg* \u\ 3 dvdv*. 
Jr n Jr n Jr n 
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By elementary changes of variables, this formula simplifies into 
with k\ = 2 p 2 N and, using f)B.3|) . 



p 2 &i / Mi,o,i(Mi i0 ,i)*|u| 3 ^d^ =2 3 / 2 p 2 &i nx 3 / 2 (Mi i o,i). 

JR n xR n 



We next observe that ^ G C°°(0,oo) and \l/ is strictly concave. It is also obvious 
that the equation ^(9) = for 9 > has a unique solution which is 9\ defined 
in fll.2ip . and that we have 

#(0) < #'(0i) - #i) = -C (0 - 0i)/2 
for some explicit constant, as well as 

(3.6) #(0) = A; 2 \el /2 - e 3 / 2 }. 



Plugging this expression for \l/ into (13.41) and using the lower bound (12. 4p on the 
temperature and the estimate (13.31) we obtain that for any e > there is C £ G 
(0, oo) such that 



(3.7) 



Voe(a 0) l) 3 / 2 - 3/2 



2+e 



<C e (l-a). 



Namely, we have thus proved that the temperature of F a converge (with rate) to 
the expected temperature 9\. In order to come back to the norm of F a — Fx, we first 
write, using Cauchy-Schwarz's inequality, 



{3.1 



l^-All^ < \\F a -M e \\ L i_ N + 1^-^^ 
< H^-MJLi+CjvllMe-Fill^, 



and we remark that 
(3.9) 



\Me-F4l, <Cp 2 3 / 2 -0 3/2 



Gathering (I3.8P with (13.91) . (13. 7p and (13.31) we deduce that for any e > there is 
C £ G (0, oo) such that 



Va G (a , 1) 



Fi|| 2 t £ <C £ p 2+£ (l-a) 



and (13.11) follows by interpolation again. 



□ 
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4 Uniqueness and continuity of the path of sta- 
tionary solutions 

Theorem 4.1 There exists a constructive ot\ G (0,1) such that the solution F a of 
U.15\) is unique for any a G [ai, 1] . We denote by F a this unique stationary solution. 

That is an immediate consequence of the following result. 

Proposition 4.2 There is a constructive constant r] G (0, 1) such that 

G, H G T a , a G (1 - 77, 1) 



implies G = H. 



Proof of Theorem 14.11 Let us assume that Proposition 14.21 holds. Then Proposi- 
tion [3J] implies that there is some explicit e G (0, 1) such that for a G (1 — £, 1] one 
has 

sup ||F Q - Fill 1,1 < V 

where r\ is defined in the statement of Proposition 14.21 Up to reducing rj, it is always 
possible to take rj < e, and the proof is completed by applying Proposition 14.21 □ 
Proof of Proposition 14.21 Let us consider any exponential weight function m with 
s G (0,1), a G (0,+oo), or with s = 1 and a G (0, 00) small enough. Let us also 
define the subvector space of L x (m _1 ) of functions with zero energy 

(4.1) = C ,o, nL 1 (m- 1 ), 

the function ip — C (\v\ 2 — N) Mi^i such that its mass is zero and its energy is 
£(ip) = 1, and LT the following projection 

niLV- 1 )^, U(g)=g-e(g)1>. 

Finally, let us introduce $ the following non-linear functional operator 

$ : [0, 1) x {W^im- 1 ) n C p>0 ) KxO, 

and 

$(1,-) : (L}(m- 1 )nC Pl0 ) -IxO, 

by setting 

$(«,</) = ({l + a)D e {g)-2N p\ IT [q q ( 9 , g) + r a A v g 

It is straightforward that &(a,F a ) = for any a G [oq, 1] and F a G T a , and that 
the equation 

$(1,0) = (0,0) 
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has a unique solution, given by g = F x = ^ denned in ffCTD . (TT2Tj) . 

The function $ is linear and quadratic in its second argument by inspection, and 
easy computations yield the following formal differential according to the second 
argument at the point (l,Fi): 



(4.2) D 2 ^(l,F 1 )h=:Ah:= [AD £ (F 1 ,h), 2Q 1 (F 1 ,h) 

where Q a is defined in ( 11. 81) and 



Ds{g,h):=bi I I gh*\u\ dvdv*. 

J JR N xR N 

Notice that we can remove the projection on the last argument in (14. 2 j) since the 
elastic collision operator always has zero energy. 

On the basis of Lemma fA.Tl in the appendix for A, we shall now prove Proposition ^. 21 
We write 

F a -H a = A~ x AF a -$(a,F a ) + $(a,H a )-AH a 

(4.3) = A' 1 (ii,I 2 ) 

with (recall that the bilinear operators Ds and Q a are symmetric) 

h := 4 D e {F u F a - H a ) - (1 + a) D(F a ) + (1 + a) D(H a ) 

h := n/ 2 ,a + n/ 2i2 

and 

^2,1 := 2 Q\(F\, F a — H a ) — Q a (F a , F a ) + Q a (H a , H a ) 
I 2i2 :=p(l-a)A(H a -F a ). 
On the one hand, 

h = 2 D(2F 1 - (F a + H a ), F a - H a ) + (1 - a) D(F a + H a , F a - H a ) 

so that 

l-^i I < C3 {\\F% — F a \\ L i + \\Fi — H a \\ L i 

+ (1 - a) \\F a \\ L i + (1 - a) ||ff a || L i) \\F a - H a \\ L x 

(4.4) < rnWWFa-HaW^-i) 

with r)i(a) — > when a — > 1 (with explicit rate, for instance 771(0;) = G\ (1 — a) 1 / 3 ) 
because of Propositions 12.11 and 13 . 11 
On the other hand, 

^2,1 = 2 Qi(-F\, F Q — iJ Q ) — 2 Q a (Fx, F a — iJ Q ) + Q a (Fi — F a , F a — H a ) + Q a (F Q — H a , F 1 — H a ) 
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From Proposition 12.11 and estimate ()A.4j) in Proposition IA.1I there holds 

Fa ~ Ha) — Qa(Fl, F a — H a )\\ L i( m -i) < C (1 — Oc) \\F a — H a \\ v^l.3( m -l) ■ 

From estimate (1A.1I) in Proposition IA.1I we have 

||<5q(-P\ — F a , F a — H a ) + Q a {F a — H a , Fi — H a )\\ L i^ m -i^ 



< C 4 [\\F a — -FiUlI^-I) + \\H a — i^lllLKm- 1 )) \\F a ~ -^allL^m- 1 )- 

Together with Propositions 13.11 we thus obtain 

(4-5) H^iIUm™- 1 ) ^ V(<X) \\ F a - #a|| L i( m -i) 

for some 77(a) — > as a — > 1. Here we can take for instance (when s = 1/2 in the 
formula of m) 77(a) = C (1 — a) 1 / 12 for some C G (0, 00) by picking a suitable e and 
interpolating. 

Finally from Proposition 12.41 there holds 

(4-6) 1 1 ^2,2 1 1 ^(m- 1 ) < C 5 (1 - a) ||F Q - I/q, I |xi( m -i). 

Gathering (14.41) . 04. 5 p and (14.61) we obtain from (14.3)) . Lemma I A. 71 and Proposi- 
tion 12.41 again 

H^a - HaW^-i) < 77(a) \\F a - HaWq^-i) 

for some function 77 such that 77(a) — > as a — ► 1 (with explicit rate). Hence 
choosing ai close enough to 1 we have 77(a) < 1/2 for any a € [ai, 1). This 

implies F a = H a and concludes the proof. □ 

Let us state now without proof, because is is completely similar to the corre- 
sponding result in [10J and is again a variation around the implicit function theorem, 
some regularity result on the path of self-similar profiles. 

Lemma 4.3 The map [ai,l] — > L 1 (m~ 1 ) ; a >— > F a is continuous on [ai,l] and 
differentiable at a = 1. More precisely, there exists F{ G L 1 (m~ 1 ) and for any 
77 G (1, 2) there exists a constructive C v G (0, 00) such that 

(4.7) \\F a - F x - (1 - a) i^|| L i (m -i) < C v (1 - a)" Va G (a , 1). 



5 Spectral study of the linearized problem 

In this section we shall enounce some results on the geometry of the spectrum of 
the linearized diffusive inelastic collision operator for a small inelasticity, as well as 
estimates on its resolvent and on the associated linear semigroup. This is straight- 
forwardly adapted from [TO] . 
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We thus consider the operator 

C*--f^Qa(f,f)+T a Af 

and some fluctuations h around the stationary solution F a : that means / = F a + 
h with / G L 1 (m _1 ) where m is a fixed smooth exponential weight function, as 
defined in (11.221) . The corresponding linearized unbounded operator C a is acting 
on L 1 (m _1 ) with domain dom(£ a ) = W^irnT 1 ) if a ^ 1 and dom(£i) = L\(rnT x ) 
(it is straightforward to check that it is closed in this space). Since the equation in 
self-similar variables preserves mass and the zero momentum, the correct spectral 
study of C a requires to restrict this operator to zero mean and centered distributions 
(which are preserved as well), that means to work in L 1 (m _1 ). When restricted to 
this space, the operator C a is denoted by C a . We denote by R(C a ) the resolvent set 
of C a , and by TZ a (i) = (C a — £) _1 its resolvent operator for any £ G R(C a ). 

5.1 The result on the spectrum and resolvent 

The result proved in this section is a translation of [TUl Theorem 5.2] for the diffusive 
inelastic Boltzmann equation. Let us define for any x G R the half-plane A^ by 

A x = {£ g C, 3?e£ > x}. 

Theorem 5.1 Let us fix ft G (fi2, 0), with /i 2 < defined in Theorem \A.6[ k, q G N 
and m a smooth weight exponential function with s G (0,1). Then there exists 
ct2 G (cti, 1) such that for any a G [ct2, 1] the following holds: 

(i) The unbounded operator C a is well defined and closed in W^'^m" 1 ) (V k > 0, 
Wq > 0). Its spectrum £(£ a ) satisfies 

where fi a is a real eigenvalue which does not depend on the choice of the space 
Wj'^m -1 ) and satisfies (with explicit bounds) 

(5.1) Ha = — 3 p (1 — a) + 0(1 — a) 2 when a — > 1. 

Moreover, \i a is a 1-dimensional eigenvalue, and more precisely, the eigenspace 
associated to fi a is R0 Q with ||0 Q ||z,i = 1 and \\<f) a — ^lll^^m-i) < C (1 — at), 
where 0i := Co (|f | 2 — iV#i) Fi t/ie "energy eigenfunction" associated to the 
linearized elastic Boltzmann operator. 

(ii) The resolvent TZ a (£) mW^m" 1 ) is holomorphic on a neighborhood of Ap\{fj, a } 
and there are explicit constants C\,Ci such that 

SUp ||fta(*)||| W *' 1 (m-i)-w5' : V-i) ^ ^ 
jzeC, 5ftez=^t 

and 

||^-a(/i + ^)|| w fe+2,i (m _ 1Hw fc,i (m -i ) < Y+UY 
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(Hi) The linear semigroup S a (t) associated to C a in Wj'^m r ) writes 



S a (t)=e^ t U a + R a (t), 

where U a is the projection on the (1- dimensional) eigenspace associated to ji a 
and where R a (t) is a semigroup which satisfies 

(5.2) II^WIIw^'^m-iHW^^m-l) < C k^ 1 

with explicit bounds. 
5.2 Sketch of proof 

The proof is straightforwardly adapted from pTJl Section 5]. We shall only mention 
the main steps of the proof and we emphasize the few points which differs here (due 
to the replacement of the anti-drift term by a diffusive term) . 

The proof is based on results on the elastic collision operator together with the 
decomposition for £ G C 

(5.3) ' £a-£ = As-B a> s 
where 

(5.4) A 5 = £+ s - £*, B a , s = v + t + (£j 3 - Cf) + P a 

with P a = L\ — C a = L\ — £a + T a A, and where C\ = C + — C* — v is the usual 
decomposition of the elastic collision operator and L\ 8 is the regularized truncation 
of the "gain" part introduced in [TTj . 

The differences with the similar decomposition in [TU] only lies in the B a § and 
P a opertors. But the key technical estimates are still true for this operator. Indeed 
it can be checked straightforwardly that the following holds (let us point out that 
one loses now two indices of regularity instead of one, which explains the minor 
change in the statement of the estimate on the resolvent). 

Lemma 5.2 (See Lemmas 5.8 & 5.9 in [10J ) Let us fix k,q > and an expo- 
nential weight function m. 

(i) There exists some constant C such that for any a G («o, 1] 

ll^«H^+ 1 2 ' 1 (m- 1 )^Vy 9 ' ; ' 1 (m-i) - ||£l _ ^a|| W |.i( m -i)_ + £i( m -i) < C (1 - a). 

(ii) There exists some constants 5* > and cti G («i, 1) such that for any £ G A w , 
5 G [0,5*] and a G [a 2 , 1], the operator 

B a ,s : < + + 1 2 ' 1 (m~ 1 ) - Wp(m- 1 ) 



is invertible and the inverse operator B a j{£) 1 satisfies 



^(m-ihWj^m-i) - dist(Ke£,Z/(R^)) 
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and 



ifrtm-iHC'K 1 ) - dist(£, i/(R*)) 
for some explicit constants C\, C% > depending on k, q, 5*, a.2- 

(Hi) As a consequence, the resolvent operator TZ a (0 satisfies for any a G [0:2,1] 
and any £ G A M2 

C 3 + C 4 II v^+ 1 1 ' 1 (m- 1 ) 

ll^«(OllwS' 1 fm- 1, i - 



im ] ~ l-C B (l-a) H^COIUm-m^-x)' 

/or some constants C{, i = 3,4, 5. 

The rest of the proof is done the same, and we recall the main steps. 
Step 1: Structure of the spectrum and rough localization 

We decompose 

l a = T + K 

with T = r Q A — L(F a ), K = 2Q+(F a ,-) - F a L(-). Since the spectrum of T 
is easily seen to be included in A^+ovi-a) (see the definition of /x 2 < in the 
statement of Theorem IA.6I and using the regularity estimates on the stationary 
solutions of Proposition 12.11) . and since the operator K is T-compact (it is the 
only point where we use the sharp estimate "in the norm of the graph" stated in 
Propositions IA.21 fc IA.3j) . we conclude (for a small enough inelasticity) thanks to 
Weyl's theorem that S(£ Q ) fl A^ only contains discrete spectrum (that is isolated 
and finite multiplicity spectrum values, or, in other words, eigenvalues), a set that 
we denote by £ d (£ a ). Moreover, thanks to point (iii) in Lemma I5T21 and to the 
geometry of the spectrum of L\ as stated in Theorem IA.6I (more particularly point 
(iii)) we deduce that £(■£«) fl A^ is confined to a disc B(0, c (1 — a)). To sum up, 
we have yet proved 

E{C a ) nA^c S d (£ a ) n B(0, c (1 - a)). 

Step 2: The "energy eigenvalue" 

First, we infer (for instance, from the decomposition (15.31) ) that any eigenfunction 
associated to an eigenvalue in A M2 is smooth: more precisely, there exists C = 
Ckqm such that V(A, with A G A M2 , ip G Wf satisfying C a ip = Xip, there holds 

U\\wZ>\ m -i)<cU\\Li- 

Second, we introduce the projector IT a on the eigenspace associated to the spec- 
trum included in A M2 . It is given by 

n a :=--L/ n a (c)dc 
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In particular the operator Hi is the projection on the energy eigenline M0i, where 
0i is the energy eigenfunction defined in Theorem IA.6I Thanks to the smoothness 
estimate on eigenf unctions just mentioned above and the estimate on C\ — C a in 
Lemma [5.21 (i) we easily deduce that 1 1 Hi — n a || < 1 for a close enough to 1, and 
then that rang IT a = rang Hi . We may sum up as 

E(£ a ) n A p = {ii a }, fi a e M, \/i a \ < c(l - a), Nul(// Q - C a ) = R(j) a . 

Third, we prove that \\<f> a — 0i|| W /M (m -i- ) < C (1 — a) exactly as in [TP] . 

Fourth and last, we obtain a first order expansion of fi a thanks to the following 
computation. By integrating the eigenvalue equation related \i a 



against \v | 2 and dividing it by (1 — a), we get 



a fa 



/; " £{<t> a ) = P 2N P { ( j )a ) - 2 (1 + a) D(F a , 



1 — a 



Using convergence of F a — * F\ and (p a — > <pi established before, and the fact that 
the mass of <f>i is zero, we deduce that 

(5.5) -^L_ E {<k) = -4 D(Fx, 00 + Oil - a). 

1 — a 



Then we compute thanks to (IB.lj) and (IB. 21) 

(5.6) 8(h) = 2N c p 6l 

where Co is still the normalizing constant in (I1.27P such that ||0i||i,i = 1. Similarly, 
using (1B.3j) . (1B.4[) and the relation fll.2ip which make a link between b\ and 9i, we 



find 

(5.7) D{FiM = lNc p 2 6 1 . 

We conclude gathering f)5.5p . (15. 6p and (15.71) . 



6 Convergence to the stationary solution 

In this section, we consider the nonlinear evolution equation (11.151) and we prove 
the convergence of its solutions to the stationary solution. 

6.1 The results 

We first state a local linearized stability result. 
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Proposition 6.1 For any a G [0K3, 1), the stationary solution F a is locally asymp- 
totically stable, with domain of stability uniform according to a G [0K3, 1). 

More precisely, let us fix p G (0, 00) and some exponential weight function m 
as in $1.22]) . There is &i,gi G N* such that for any Mq G (0, 00) there exists 
C,eG (0, 00) such that for any a G [03, 1], for any f m G H kl D L 1 (m~ qi ) with mass 
p, momentum satisfying 

(6-1) ll/in|| J H-fcinL 1 (m~' Jl ) < M), II /in ~ ^alU^m" 1 ) < £ > 

i/ie solution f to the equation $1.150 with initial datum f- in satisfies 

(6.2) Vt > 0, ||n Q (ft - F )|| L i (TO -i) < C ||/in - Fall^i^-i) e**"', 

(6.3) Vt > 0, || (Id - n a ) (/ t - F a )|| il(m -i } < C ||/ in - F a || £ i (m -i) e ( 3 /2)^*. 

Then we prove that when the inelasticity is small, depending on the size of 
the initial datum (but not necessarily close to the stationary solution), the equa- 
tion (11.151) is stable. This mainly relies on the fact that the entropy production 
timescale is of a different order (much faster) that the energy dissipation timescale 
as a — > 1. 

Proposition 6.2 Define k 2 '■— max{fc ,^i} ; Q2 '■= m ax{g , <Zi, 3}, where ki and 
q, t are defined in Theorem \A.5\ and Corollary \A.4\ For any p, 8 , M there exists 
«4 G [oj 3 , 1), Ci G (0, 00) and for any a G [04, 1] there exist <p = <p(a) with cp(a) — > 
as a — > 1 and T = T(a) (possibly blowing-up as a — > 1) such that for any initial 
datum < / in G L\ 2 n H k2 n C p>0) e with 

ll/in|Ll 2 nH fe 2 < M , 

the solution f associated to the rescaled equation $1.24 ) satisfies 

vt > 0, s(f t ) > a 

and for all a' G [0:4, 1) and then all a G [a', 1] 

(6.4) Vt>T(o/), \\ft-F a \\ Ll <¥>(«')• 

11 hjj 2 

Then the proof of the gobal convergence for smooth initial data only amounts 
to connect the two previous results of Propositions 16.11 and 16.21 by choosing a such 
that (p(ct) < e where e is the size of the attraction domain in Proposition 16.11 and 
(p(at) is defined in Propositions 16.21 More precisely, we have straightforwardly the 

Corollary 6.3 Let us fix an exponential weight function m as in $1.22\) . with ex- 
ponent s G (0,1). Then for any p, £ Q , M there exists C and o 5 G [0:4,1) (de- 
pending on p, £ Q , M ,m) such that for any a G [05,1) and any initial datum 
< /n G L 1 {m~ q2 ) H H k2 satisfying 

/in G C p fl,£ , ll/in||L 1 (m-'J2)n// fc 2 < M , 



23 



the solution f associated to the reseated equation {1-2$ satisfies 

Vt>0, ||n a (/ t -F a )|| L i (m -i) <Ce^\ 

Vt > 0, ||(Id - n Q ) {f t - F a )|Ui (m -i) < Ce^^K 

As a by-product of the previous propositions, we state and prove a result which 
provides a partial answer to the question (important from the physical viewpoint) 
of finding Liapunov functionals for this particles system. Let us define the required 
objects. We consider a fixed mass p and some restitution coefficient a whose range 
will be specified below. At initial times, non-linear effects dominate and therefore 
we define 

«!(/) :=H(g\M[f]) + (S-S a ) 2 

where £ a = S(F a ) is the energy of the self-similar profile corresponding to a and 
the mass p. At eventual times, linearized effects dominate. Therefore we define a 
quite natural candidate from the spectral study: 

f + CO 

: = OTlW) + (!-«)/ h% ds, 

Jo 

with h 1 = H a h, h 2 = H^h and h = g — G a . 

Proposition 6.4 There is k^ G N big enough (this value is specified in the proof) 
such that for any exponential weight function m as defined in U.2S\) . any time 
t G (0, oo) and any p,S ,M G (0, oo), there exists G (0, oo) and a 6 G [a 5 , 1) 
such that for any a G [a^, 1] and initial datum f in G H ki H L 1 (m _1 ) satisfying 

fin e C p>0 , So , \\fin\\ H HnL Hm -l) < M , fin(v) > M^ 1 e~ M ^\ 

the solution g to the rescaled equation with initial datum f in is such that the 

functional 

H(f t ) = Hi(ft) l{ Wl(/t) > K ,} + H 2 (ft) 1 { Hl (f t )<n,} 

is decreasing for all times t G [0, +oo). Moreover, H(f(t, •)) is strictly decreasing as 
long as fit,-) has not reached the self-similar state F a . 

Then the previous results can be extended to general initial data by a study of 
the decay of singularities in the same spirit as in [10]. Indeed one can easily prove 
the following decomposition result. 

Lemma 6.5 Consider f in G L\ and the associated solution f G C([0, oo); L|) to the 
rescaled equation \1.21$ . Assume that for some constant p,c\,M\,T G (0, oo) there 
holds 

(6.5) / m eC p , , H/inlUi < M l5 ViG[0,T], £(/(*,-))> ci. 

Then, there are a 7 G [a 6 , 1) and X G (— oo, 0), and for any exponential weight 
function m (as defined in M.22\) and any k G N, there exists a constant K ( which 
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depends on p, c\, M 1; k, m) such that for any a G [a 7 , 1], we may split f = f s + f R 
with 

(6.6) Vt G [0,T], ||/^,.) |U fc nLHm-i) < K, \\f R {t,-)hl < K e xt . 

This last tool together with the L 1 usual stability result of the elastic Boltzmann 
(stating that the error between two flows grows at most exponentially) allow to 
conclude to the proof of point (iv) of Theorem 11.31 

6.2 Sketch of the proof 

The proof is exactly the same as in [TUl Section 6]. Indeed the precise form of the 
equation was not used in this proof, it only uses: 

• on the one hand, for the perturbative argument, the spectral properties ob- 
tained before (and valid here), some uniform estimates on moments and reg- 
ularity in terms of a lower bound on the energy (still valid here) and some 
estimates from above on the bilinear term (which are the same here); 

• on the other hand, for the argument "in the large", it uses elastic entropy - 
entropy production estimates (independent of our problem), and the difference 
of timescales between entropy production (which is 0(1)) and temperature 
thermalization (which is 0(1 — a)), also present here. 

A Appendix: Functional toolbox on the collision 
operator 

In this first appendix we shall recall some functional results on the collision operator 
Q which were obtained in [TO] . 

First let us show that the collision operator depends continuously on the in- 
elasticity coefficient a G [0,1]. Since it is an unbounded operator, this continuous 
dependency is expressed in the norm of the graph of the operator or in some weaker 
norm. We start showing that this dependency of the collision operator is Lipschitz, 
and even C l ' n for any 77 G (0, 1), when allowing a loss (in terms of derivatives and 
weight) in the norm it is expressed. Let us define the formal derivative of the collision 
operator according to a by 

Q' a (gJ):=V v - ( [ [ g(%(a))f(>v(a))b\u\ (^P) dadv* 
or by duality 

(Q'a(9,f),il>) : = [ [ [ 9*fb\u\ (^-^ — -J VTp(v' a )dadv*dv. 
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Proposition A.l (See Proposition 3.1 in |10| ) Let us fix a smooth exponential 
weight m = exp(— a \v\ s ), a G (0, +oo) ; s G (0, 1). Then 

(i) For any k, q G N the exists C G (0, oo) such that for any smooth functions f, g 
(say in S(M. N ) ) and any a G [0, 1] there holds 

(A- 1 ) HG«(0>/)llw*' 1 (m-i) ^ C *>rn \\f\\w^\(m^) IMIw£i(n»-i) 

( A -2) \\Q' a (9, f)\\ W ^(m-^ < G k,m \\f\\w*+l'\m-i) IMI W^'V*- 1 )" 

(ii) Moreover, for any smooth functions f, g and for any a, a' G [0, 1], there holds 

\\Qt(g, f) - Q$(jg, /)-(«- «') c&fo, /) ll^v- 1 ) 

( A -3) < | a - o/\ 2 ||/IUj +8 ( m -i) lblUi +3 (m-i)- 

fwj ^4s a consequence, there holds 

\\Qt>{g, f) - Ql(9, /)llw*( m -i) <C\a- a'\ WgWw^im-iy 
(A.4) 

and for any 77 G (1, 2), there exists k v G N ; q v G N and G (0, 00) sitc/i £aai 

WQUg, f) - Qa>(g> /)-(«- «') tfafo /) 11^-) 

(A.5) < C, |o; - a/|" H/H^.i^ lbll^ (m -i)- 

We next state a mere (Holder) continuity dependency on a, which is however 
stronger than Proposition IA.1I in some sense, since it is written in the norm of the 
graph of the operator for one of the arguments. 

Proposition A. 2 (See Proposition 3.2 in |10j ) For any a, a' G (0,1], and any 
g G L\(m~ l ), f G W^ 1 (m~ 1 ), there holds 

[ \\Qt{g, f ) - Qt'(g> /)ILi(m-i) ^ £ ( a ~ a ') 11/11 W^fm-i) lblUi(m-i)» 

(A.6) { 

[ \\Qt(f,g) - <5a'(/^)llLi(m-i) < £ (« - «') ll/IU^Cm-i) H# II ^(m" 1 ) " 
i 

where e(r) = C r 3+4 / s /or some constant C (depending only on b). 

As a consequence of Propositions IA.1I & IA.2| together with Lemma 14.31 and 
Proposition 12. 1[ we have the 

Proposition A. 3 (See Proposition 5.8 in |10j ) For any k, q > and for any 

exponential weight function m there exists some constant C such that Va G (ojq, 1] 



(0 1 








(ii) 1 




- £+1 




< C(l-a 


m) 








< £(1 - a), 



where e is defined as previously (up to a constant). 
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Second we shall give some estimates on the entropy production functional asso- 
ciated with the elastic collision operator Qi. We begin with a simple consequence 
of Proposition IA.1I 

Corollary A. 4 (See Corollary 3.4 in |10j ) There exists k ,q G N such that for 
any a,i G (0, oo) % = 1, 2, 3 ; there exists an explicit constant C G (0, oo) such that 
for any function g satisfying 

\\g\\ H ^L\ o < ai, g> a 2 e- a3lvl \ 

there holds 

\D Hja {g)-D Htl {g)\<C{l-a), 
where we recall that D H a is defined by 

(A.7) D H>a {g) = l Iff gg (lA-i og lA-i\ Bdvdv*da>0. 

2 J J J R 2n xS n-i \gg* gg* J 

Let us recall now two famous inequalities, namely the Csiszar-Kullback-Pinsker 
inequality (see [HIE]) and the so-called entropy-entropy production inequalities (the 
version we present here is established in [13]) that we will use several time in the 
sequel. 

Theorem A. 5 (See Theorem 3.5 in [10J ) (i) For a given function g G L\, let 
us denote by M[g] the Maxwellian function with the same mass, momentum 
and temperature as g. For any < g G Ll(W N ), there holds 

(A.8) \\g - M[g]\\ 2 Ll < 2 p(g) [ gfo-i-dv. 

J RN M [g\ 

(ii) For any e > there exists k £ , q £ G N and for any A G (0, oo) there exists 
C £ = C £ ^a £ (0, oo) such that for any g G H ke n Ll e such that 

g(v) > A' 1 e~ A ^\ \\g\\ HhsnLl < A, 

there holds 

(a.9) c e , (J) .-.^ Jln _|_ (it ,y +e <z, H , 1 ( 9 ). 

Third and last we recall some results on the linearized elastic collision operator. 

Theorem A. 6 (See and Proposition 5.5 & 5.7 in |10j ) Let define C\ as 
the restriction of C\ to the space L\ n C 0i o- Then for any k, q > and for any 
exponential weight function m the following holds. 
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(i) The regularized truncation Cf s of the "gain" part (see (f5.^| ) and fTDj) satisfies 

~ £ MlU^( m -l HW /M (m -l) ^ £ ( 5 ) 

where e(5) > is an explicit constant going to as 5 goes to 0. For any 5 > 0, 
the linear operator As (see 115. 3\) ) satisfies As : L 1 — > W^' l {m~ l ) is bounded. 

(ii) The unbounded operator Ci is well defined and close in Wj'^m" 1 ). Its spec- 
trum is real and satisfies 

S(A) C A% U {/ij, 

with /i2 < and {i± = is the 1- dimensional "energy eigenvalue" associated to 
the "energy eig en function" 0i := Co (\v | 2 — iV^i) iq. In particular, L\ is onto 
from O fl Lj(m _1 ) onto (9. 

(Hi) The resolvent TZ\{^) has a sectorial property for the spectrum sub str acted from 
the u energy" eigenvalue, namely there is a constructive A G (/U2,0) such that 

V£ G A H^OIIw^W) < «M + 



some explicit constant ak, q , bk, q > and 

37T 37T 



^ = U 6 C, arg(£ + A) G 



4 ' 4 



and < — 



A quite simple consequence of point (ii) above is the following quantitative in- 
vertibility result on a modified version of the linearized elastic collision operator. 

Lemma A. 7 (See Lemma 4.3 in [10J ) The linear operator A : Lj(m _1 ) — > 
1x0 defined thanks to \4-W ^ s invertible: it is bijective with A^ 1 bounded with 
explicit estimate. 



B Appendix: Moments of Gaussians 

We state here some results on the moments of tensor product of Gaussians (the 
proof is done is the appendix of [TO]). 

Lemma B.l The following identities hold 
(B.l) / M 1)0 ,i \v\ 2 dv = N, 

(B.2) / M 1A iM 4 dn = A (A + 2), 
Jr n 

(B.3) / M lfitl (M lfiA )*\u\ 3 dvdv* = 2 3/2 [ M lfiA \v\ 3 dv, 

JR N xR N JR n 

(B.4) f Afi, 0l i (Mi )0 ,i), \v\ 2 \u\ 3 dv dv* = V2(2N + 3) [ M lfi ,i{v) \v\ 3 dv. 

JR N xR N JR n 
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C Appendix: Interpolation inequalities 



Again the proof of the following simple interpolation inequality can be found in the 
appendix of [TO] . 

Lemma C.l (i) For any k, k*,q,q* G Z with k > k* , q > q* and any 9 G (0, 1) 
there is C G (0, oo) such that for h G W q »* ^{mT 1 ) 

fc**,g** G Z such that k = (1 - 6) k* + 6k**, q = (1 - 6) q* + 6q**. 

(ii) For any k,q G iV* and any exponential weight function m as defined in M.22fy . 
there exists C G (0, oo) such that for any h G H k% fl L 1 (m -12 ) with : = 
8fc + 7(1 + N/2) 

D Appendix: Lower bound for the diffusive in- 
elastic Boltzmann equation 

In this last appendix, we state and prove a technical result on the pointwise lower 
bound of the diffusive Boltzmann equation (the new difficulty as compared to pre- 
vious results is due to the diffusion and not transport nature of the term added to 
the collision operator). 

Proposition D.l Let g G C([0, oo); L\) be a solution of the rescaled evolution equa- 
tion 

d t g = Q(g,g) + r a A v g, 

with inelasticity parameter a G (0, 1), and assume that for some C,T G (0, oo) 

sup||#|| L 2 nL i < C. 

[0,T] 

For any t\ G (0, T) there exists a\ G (0, oo) (depending on C, p and t\ but not on 
T) such that 

(D.l) Vt G [h,T], Vv G R N , g(t,v) > a^ 1 e~ a ^\ 

We closely follow the proof of the Maxwellian lower bound for the solutions of the 
elastic Boltzmann equation (see 0,[12]) taking advantage of some technical results 
established in its extension to the solutions of the inelastic Boltzmann equation (see 
[9j Theorem 4.9] and pJJ Lemma 2.6]). The starting point is again the evolution 
equation satisfied by g written in the form 

(D.2) d t g - r a A v g + X(v) g = Q+(g, g) + (X(v) - L(g)) g, X(v) = k(1 + \v\), 
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where the last term in the right hand side term is clearly non-negative for some well- 
chosen numerical constant k G (0, oo). Let us introduce the semigroup S t associated 
to the linear evolution equation d t g — r a A v g + X(v) g = 0. 

We start establishing some technical results that we need in the proof of Proposi- 
tion [DJ3 Let us recall some elementary results extracted from [9] (which are a mere 
adaptation to the inelastic collision of some result proved in [12]). 

Lemma D.2 J3 Lemma 4.6] Let < £ G {L\ n L 2 ){R N ) satisfy 
£dv = mo, / £\v\ 2 dv<mi, / £ 2 dv < 



for some positive real constants mj. There exists R > r > and 77 > depending 
only on m ,m 1 ,m 2; and (vi)i=i,...,4 such that \vi\ < R, i — 1, . . . , 4, \v j — Vj\ > 3r 
for 1 < i 7^ j < 3, and 



(D.3) / £(v)dv>7] for i = 1, 2, 3, 

JB( Vi ,r) 

(D.4) Vwj G -B(fj, r), E% Wi fl W2 is a sphere of radius larger than r, 

where, as in [9, Proposition 1.5], E° 3W4 stands for the hyperplan orthogonal to the 
vector W3 — uii and passing through the point fl(w 3 , W4), defined by 

0(w 3 ,u7 4 ) := w 3 + (l- a) (w 3 - w 4 )/(l + a), 

and W2 stands for the sphere of all possibles post-collisional velocity v' defined by 
(O) from (v = w 3 , v* = u> 4 , a) . 



Lemma D.3 Let < f,g,h G L 1 (M N ) and (t>i) i= i,...,4 satisfy / POP . (LT$) and 



\vi\ < R, i = 1, . . . , 4, \vi — Vj\ > 3r for 1 < % ^ j < 3 ; for some given constants 
R > r > 0, r] > 0. There exists Tq > 0, S' Q > and r)' > only depending on R, r 
and 7] such that 

Vt'e [0, T ] Q+(f, S t ,Q + (g, h))> V ' 

Proof of Lemma \D.3\ We first recall a convenient formula to handle representations 
of the iterated gain term. Using Carleman representation [91 Proposition 1.5], for 
any /, h and i and any v G M. N there holds (setting 'v = w and 't>* = w*) 

Q + (f,S t Q + (hJ))(v)=C' b [ -r^-rl f S t Q + (g,h)(w*)dwA dw 

J R N \V-W\ [J Evw J 

f(w)g(z)h(z*) \ C' h ^ 4 / ( s tlEzJ(z')b((T ■ z)dcr\ dz^dzdw, 

{ \V-W\ J s n-i v ' w J 

where E% w is defined in the statement of Lemma ID. 21 and z' is defined from (jl.4p 
with (v = z,v* = z*,cr). Let us define / := f l B (v 4 ,r), 9 = glB( Vl ,r), h = hl B (v 2 , r ) 
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and S t the semigroup associated to evolution equation d t g — r a A v g — \(2R)g = 0. 
By the maximum principle, we have for any v G B(v^, r) 

Q + (f,S t Q + (h,£))(v) > [ f(w)g(z)h(z*) [c'^- [ (S tX )(z')da) dz*dzdw. 

J R 3N I ti J S N-1 J 

Taking v G B(vs, r), w G B(v±, r), z G B(vi, r), G B(v2, r) and denoting by A the 
term between brackets, we have thanks to flD.4j) 

A{v,w,z,z*):=CZ?- [ [ X (z' ' -u)e-^dadu>^-C b Cr N - 2 l, 
R Js N ~ 1 Jr n k * 

for any a G [0,1] and t G [0, T ] with T small enough. We then conclude the proof. 
□ 

For a given fi G (0, 1) we consider a function \ G W /2,00 (M + ), which satisfies the 
following properties: x > 0, x' < 0, x = 1 on [0, yu] and x(a;) = (1 — x) 2 for any 
x G [//, 1] with // := max{l — 1/iV, (l+/i)/2}. Abusing notations we define the radial 
symmetric functions x an d X s on ^ N by setting x(t>) = an d — 

for any 5 > 0. 

Lemma D.4 For any 5o > 0, there exists C G (0, oo) such that 

V«G(0,1), Vi?>0, VveB(0,R), V5>5 , Vt > 0, 5* fox,) > e " c ( 1 + fl + a )* ( % ). 

Proof of Lemma ID. 41 We shall rely on a maximum principle argument for the 
operator 

Tf:=T a A v f-X(v)f. 

First it can be seen easily that if fo is a non-negative function and / = f(t,v) 
satisfies the differential inequality 

d t f>Tf, tG[0,T), 

then ft remains non- negative on this time interval [0, T). 
Second let us show that 

(j)(t,v) := e~ at TvXs 

is a sub-solution for a = as^R big enough, on a small initial time interval [0, T), in 
the sense that 

d t <p<T<p, te[0,T). 

After elementary computations and taking advantage of the fact that x is radially 
symmetric, it amounts to show that 

(D.5) [a-X(5w + v)} X +^ 



ft i 

X + 



(N-l) 



\w\ 



X 



> 



on 
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for a big enough. From the definition of X , it is zero outside of -8(0, 1) and therefore 
equation (ID. 51) holds on M. N \B(0, 1). Observing that 

x "( x ) + i^l x '( x )>0 
x 

for any x G [//, 1], we deduce that equation (ID. 5[) holds on B(0, 1)\B(0, //) as soon 
as a > k (1 + 5 + R). Finally, using that X is decreasing and bounded in W 2 '°°(0, 1), 
we see that equation (ID.5j) holds on B(0, fjf) as soon as 



[a-K(l + 6 + R)] X (jS) > 



5 2 



x »( x ) + VLJl x >( x) 



X 



L°°(0,1) 



and therefore as soon as a > C (1 + S + R) for some constant C G (0, oo) only 
depending on X , <5 , p and N. □ 

Proof of Proposition ID. 11 We split the proof into two steps. 

Step 1. Thanks to the Duhamel formula a solution g to the evolution equation (1D.2j) 
satisfies t 

(D.6) g(t,-)>S t g(0,-)+ [ S t ^ s Q + (g(s, ■), g(s, ■)) ds. 



Let us fix t > and define g (A ■) := 9 (to + •)• Using twice the inequality (ID.6I) . 
we find 



(D.7) g (t,-) > 



o Jo 



S t ^ s Q + (S s g , S s ^ s iQ + (S s ,g , S s >g )) ds' ds. 



We first apply Lemma ID.3I to i = S s > go with r = s — s' and we next use Lemma ID.4I 
to obtain for any t G [0, T ], with T given by Lemma ID. 31 



t rs 



(D.8) g (t,-) > 




o Jo 



S t - s r}l B (v,r)ds' ds>rj [ e a(1 (t s)) r dXr sds. 

Jo 



We have then proved that there exists 7\ > and v G -8(0, R) and for any t\ G 
(0, Ti/2] there exists rji > such that 

Vt G [0,Ti/2], <h(f,0 := <7o(* + *i,0 > Vi T vX Sl =■ Vi fa- 
Step 2. Using again the Duhamel formula (1D.6|) and the preceding step we have 

9i(t,-)> / S t _ s Q + (g 1 (s,-),g 1 (s,-))ds. 
Jo 

Now let recall, that on the one hand, from [21 Lemma 4.8], there exists k G (0, oo) 
such that 

Qa(l-B(0,l)) 1-B(0,1)) > «' lfi(0,V5/2) 
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and that, on the other hand, the scaling properties of the term infers that 

QZMM *(./*))(*) = ^ QM 

for any function <p and scaling coefficient 5 > 0. Then, thanks to Lemma [D .41 there 
holds 

01 (V) > »7i J S t - s Q + (x h ,X Sl )ds 

on [0, T 2 ] with T 2 G (0, Ti/2], 9 G (1, y/E/2) as close as we wish to a/5/2 by choosing fi 
close to 1 and b a numerical cosntant (depending on k', x)- Defining t k = t^-i+^i 2~ fc 
and repeating the precedent computation we see that 

g{t,-)>VkX Sk Vte[t k ,T], 

with 5k = k 5\ and 

Vk+x = vlbe-^+^^^^S^ 1 (6- N ^) k = 7]l A B k , 
for some constants A,B G (0, oo). Elementary computations yields to 

for some constant D G (0, oo). As a conclusion, using that 6 8 > 2, we have proved 

Vt > ti, Vfc G N, #(t,<u) > D 2 * l B(fi)2fc /8 (<5l/2)) (w), 

from which we easily conclude. □ 
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